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A phenomenological model for the quark structure of 
mesons is considered. The model is based on the tube model 
for QCD, where all quanta with nonzero transverse momenta 
are neglected. In the limit that the mass term of the gluons 
goes to infinity, the model is equivalent to a combination of 
the 't Hooft and Gross-Neveu models and can be solved semi- 
analytically. The model has the properties of confinement, 
chiral symmetry breaking and asymptotic freedom and thus 
resembles QCD in three key respects. Spectra, distribution 
amplitudes and form factors of mesons are analyzed. 



I. INTRODUCTION 

Many high energy scattering processes, probe the 
quark-gluon structure of hadrons predominantly in one 
direction It is thus suggestive to consider a phe- 
nomenological model for QCD, where the transverse mo- 
menta of all constituents are neglected. At a semi- 
classical level, one can think of this approximation as a 
formulation of QCD in a "tube" with periodic boundary 
conditions. By taking the radius of the tube very small, 
all modes other than /cj^ = have a very high energy 
and thus freeze out. The model which will be considered 
in this work is based on such a classical reduction of the 
QCD Lagrangian to an effective 1-t-l dimensional theory, 
which is subsequently quantized. It should be empha- 
sized that the model thus obtained is not a dimensional 
reduction at the level of the full quantum theory, which is 
slightly more subtle. As a consequence, the model should 
not be considered a rigorous approximation to full QCD 
but rather a phenomenological model. Furthermore, it 
should be noted that even though the model is formu- 
lated in 1+1 dimensions, one should not consider it a 
mere toy model but rather a simple (because all fcj^ van- 
ish) phenomenological model for QC-D3+1. 

While it may be possible to understand many features 
of high energy scattering experiments in terms of a sim- 
plified effective theory, which is reduced to degrees of 
freedom with vanishing transverse momenta, boost in- 
variance often seems essential for microscopic descrip- 
tions of such experiments. Furthermore, high energy 
scattering experiments typically probe correlations along 
light-like directions 0|. It thus makes sense to employ 
light-front (LF) quantization (2|-|^] when studying dimen- 
sionally reduced QCD. 

Such a model has been considered in Rcf. Q for glue- 
balls and in Ref. M for mesons, where discrete light-cone 



quantization (DLCQ) Q has been used as a numerical 
tool. However, due to the rather singular behavior of 
the fermion wave functions, the resulting spectra are not 
very well converged numerically when the quark masses 
become small. For form factors, which we would like to 
consider in this work, it is even more important to de- 
scribe the end-point behavior of the wave functions very 
accurately. Improvement methods based on DLCQ ||] 
have so far only been developed for gluons. 

In order to develop a model which can be solved semi- 
analytically, we will consider a modified version of the 
model in Ref. |Q where we assume that the mass term for 
the effective gluon is very large and thus gluon degrees of 
freedom can be eliminated perturbatively. The resulting 
effective interaction for the quarks resembles 't Hooft 's 
large- iVc QCD in 1+1 dimensions |jl^ with an additional 
helicity dependent Gross-Neveu interaction |jll|,|l2| 

Because the matter degrees of freedom in collinear 
QCD couple to a 1+1 dimensional gauge theory, con- 
finement is an almost trivial feature of the model. The 
effective Gross-Neveu interaction, which arises from elim- 
inating the transverse gluons, provides an induced mass 
for the quarks and thus leads to spontaneous breaking 
of chiral symmetry. Furthermore, the Gross-Neveu in- 
teraction is also asymptotically free. Thus, even though 
collinear QCD cannot be derived rigorously as an ap- 
proximation to full QCD, the mere fact that it shares 
with real QCD the important properties of confinement, 
chiral symmetry breaking as well as asymptotic freedom 
makes it worth while to investigate collinear QCD as a 
phenomenological model. 

As a first application of the model, we will study the 
low lying meson spectrum. This part of the calculation 
will also be used to determine the free parameters of the 
model. Once all terms in the Hamiltonian are fixed, we 
will then proceed to calculate other observables, specifi- 
cally distribution amplitudes (i.e. light-cone wave func- 
tions) and form factors. 



II. THE COLLINEAR QCD MODEL FOR 
MESONS 

The basic idea of collinear QCD (also referred to as 
the tube model) is to start from the classical QCD La- 
grangian and to neglect all dependences on transverse co- 
ordinates (transverse with respect to an arbitrarily cho- 
sen, but fixed direction), yielding 
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where we have introduced LF coordinates A^ = A^ — 
ztA^y To keep the calculations simple, which will 
help in identifying the essential physics, only the limit 
Nc 00 will be considered here. 

Neglecting modes with k± breaks invariance un- 
der x± dependent gauge transformations. Therefore a 
mass for the transverse component of the gauge field has 
been added, sinc e it is now no longer protected by gauge 
invariance. Eq. (2.1) is still invariant under gauge trans- 
formations that depend on only and we are thus free 
to choose the gauge — 0. [|. Neither A^ nor 
(where ip± = ^7^7^V') &re dynamical degrees of free- 
dom, since their LF-time (x"*") derivative does not enter 
the Lagrangian (2.1). This degrees of freedom are elim- 
inated from the Lagrangian [Eq. (^.l|)] using the cor- 
responding constraint equations (Euler Lagrange equa- 
tions) before quantizing the theory. 

The dynamical degrees of freedom of the model are 'ip+ 
and A±, which are quantized canonically. In summary, 
what the dimensional reduction yields is a 1-1-1 dimen- 
sional gauge theory, coupled to both scalar matter in the 
adjoint representation {A±) and fermionic matter in the 
fundamental representation {ip). Both A± and V' have 
internal helicity degrees of freedom — remnants of the 
underlying 3-f 1 dimensional theory. 

In addition to the coupling of the matter fields to the 
1-1-1 dimensional gauge theory, there is also a direct, 
Yukawa-like, coupling between fermions and A±^ which 
flips the helicity of the fermions. 

Note that the dimensional reduction procedure, the 
choice of gauge (A+ = 0) as well as the quantization 
plane (x~^ — 0) are all manifestly invariant under ro- 
tations about the z-axis. As a result, the z-component 
of the total angular momentum is a conserved quantity 
in the model. Since there is no orbital angular momen- 
tum left after the dimensional reduction, this means that 
the total z-component of the spin, i.e. the sums of z- 
components of the spins of quarks, plus anti-quarks plus 
gluons, is conserved. This feature of the model will be 
very helpful when it comes to classifying states. 

Such a model has been studied in Ref. 0], where ap- 
proximate numerical solutions have been obtained using 
a Fock space expansion and DLCQ ||^. In Ref. , one 
can also find an explicit expression for the dimensionally 
reduced LF Hamiltonian for QCD, which will not be re- 
produced here in its full generality, since we will consider 
only a solvable limiting case of the model in this work. 



■0 By allowing the mass term for A± to become arbi- 
trarily large (with appropriate rescaling of the coupling 

2^ constant in order to keep the physics nontrivial), it is 
possible to systematically eliminate A±^ yielding an ef- 
fective interaction which acts only on the fermion degrees 
of freedom. For this we note that, as long as they are not 
diverging, all interactions of Aj_ become negligible in this 
limit. Furthermore, we can neglect all interactions in in- 
termediate states which contain quanta of A±_ and are 
thus highly off-shell. This makes it an easy task to elimi- 
nate A± using a LF Tamm Dancoff procedure. First this 
gives rise to a kinetic mass counter-term for the fermions. 
As is explained in the Appendix, the Tamm-Dancoff pro- 
cedure also requires to renormalize the one giuon vertex. 
After eliminating all gluon dressing from the quark lines, 
we are thus left with a fermion field that has a mass term 

M. g 

In order to determine the effective qq interaction aris- 
ing from eliminating A±, let us consider qq scattering 
(Fig. 0)0 




x-y 




FIG. 1. Scattering of a quark with initial LF-momentum x 
from an anti-quark of initial momentum 1 ~ x. The dotted 
line represents the exchanged gluon. Only one time ordering 
is shown. 

Since the quark helicity flips at the quark gluon vertices 
and since total (quark+gluon) helicity is conserved at 
each vertex, the one gluon exchange interaction vanishes 
identically for helicity ±1 mesons, i.e. when quark and 
antiquark have the same helicity. Thus, as our first result, 
we find that mesons with helicity ±1 are in our model 
described by 't Hooft's equation only 



^Explicit zero-modes, i.e. modes that do not depend on x 
will be neglected throughout this work. 



Strictly speaking we would have to keep separate kinetic 
and vertex mass terms. However, as is discussed in Appendix 
H, choosing them equal gives rise to a consistent solution. 

Once again, since A 00, all other interactions can be 
neglected in intermediate states with A± quanta. 
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The only role played by the transverse gluons for these 
mesons is to give the quarks a "constituent mass" M . 

The situation is different for mesons with helicity zero, 
i.e. when quark and antiquark have opposite helicity, 
in which case the one gluon exchange matrix element is 
nonzero and reads for the time ordering depicted in Fig. 

m 
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The other time-ordering gives (for A — > cxd) the same 
expression, but with 0(y — x), so that the 0-function 
disappears in the sum of the two time-orderings 
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where the different Lorentz structures appearing in Eq. 
(2.3) have been separated. Let us first consider the 
"pion" channel, i.e. mesons with helicity zero, an odd 
spin wave function and an even orbital wave function 
for the quarks. In this channel, the effective interaction 
reads 
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(2.5) 



where 5e// ~ 9^ 1^ The first interaction term in Eq. 
(2.5) is the same as what one gets in the chiral Gross- 
Neveu model, i.e. an attractive s-channel ipii"f5ipi'>Pji"f5tpj 
interaction, where i,j are color indices. Such an interac- 
tion leads to a divergence when it is iterated and one has 
to introduce a cutoff. From studies of the Gross-Neveu 
model , it is known that the the coupling g^f / goes 

to zero as the cutoff is sent to infinit y (as ymptotic free- 
dom!). The other two terms in Eq. (2.5) are what one 



gets if one introduces an s-channel vector interaction of 
the form ilii'^^tliiipj^^i^j in 1 -|- 1 dimensions. Such an in- 
teraction does not lead to divergences (in covariant per- 
turbation theory) — neither when iterated with itself nor 
in conjunction with the QCZ^i+i or Gross-Neveu inter- 
actions. Thus as the running coupling goes to zero, this 
s-channel vector interaction becomes less and less impor- 
tant, until it disappears when the cutoff goes to infinity. 
In the effective interaction, one can thus omit this term 
entirely, yielding 
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(2.3' 



In the scalar channel (total helicity zero, spin part odd. 
Orbital part odd) one finds the same result as for the 
;hiral Gross-Neveu model in the scalar channel 
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(2.7) 



and also the relative sign and strength between the tt 
and a channel are consistent with the chiral Gross-Neveu 
model. 

In addition to this contact interaction from eliminating 
the transverse Gluon field, the mesons are also governed 
by the QCDi+i interactions. Thus, after all the dust has 
settled, we are left with an effective two dimensional field 
theory, describing mesons both in the tt and a channels, 
which in covariant notation reads ^ 



C 
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(2.8) 



where ^/^ is a two-spinor with N color components 
(fundamental representation), F^^, = d^A^, — d„Af^ + 
i^L [Afj^,A^] is the gauge field and the hmit iV — + oo 
is implied. Obviously, for 7 = one obtains the 't Hooft 
model {QCDi+i in the hmit iV ^ 00) (lO), 

Finally, let us turn to mesons with zero total helicity, 
but even helicity wave functions. By the same reasoning 
as above, one obtains the same interaction as in the n — a 
channels, but with the sign reversed (because the inter- 
action is a spin exchange interaction), i.e. they act repul- 
sively. In non-perturbative calculations, there is a crucial 
difference between renormalized zero-range interactions 
that are attractive and those that are repulsive! This fact 
is best illustrated in the context of a simple example, the 



* The overall coupling is irrelevant at this point, since it is 
renormalized anyways. 



^One can arrive at the same result also by means of a rather 
lengthy Fierz rearrangement. However, it is much more illu- 
minating to derive this result in the LF helicity basis as has 
been done above. 
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non-relativistic Schrodinger equation in two spatial di- 
mensions with a delta function interaction. Since a delta 
function in two dimensions leads to a logarithmic diver- 
gence, one needs to introduce a cutoff. Suppose now that 
the delta function is attractive. Then, in order to keep 
a finite mass of the bound state, the coupling constant 
must go to zero as the cutoff goes to infinity — otherwise 
the wave-function gets sucked into the singularity. Now 
suppose one calculates the effect of a repulsive interac- 
tion with the same strength but opposite (repulsive) sign. 
In this case, there is no enhancement of the interaction 
due to being sucked in and as a consequence the repul- 
sive effects of the interaction completely disappear as the 
cutoff is sent to infinity. 

The same happens here with the spin-even helicity zero 
interaction. Because of the symmetries of the underly- 
ing action, the coupling constant is the same as in the 
attractive case (tt — cr channel) but of opposite (repul- 
sive) sign. Since non-perturbative renormalization for the 
TT — cr channel (as done in Appendix ^) yields a vanish- 
ing bare coupling, and since the same bare coupling acts 
in the repulsive case, one finds that the effective interac- 
tion from one gluon exchange in the repulsive channels 
completely disappears. The spin-even helicity mesons 
are thus described by 't Hooft's equation (without Gross- 
Neveu interaction). This is a very interesting result, since 
it implies that the spin-even helicity mesons are de- 
scribed by the same bound-state equation as the helicity 
±1 mesons, and they thus form degenerate triplets — 
just as one would expect from vector and axial vector 
mesons! 

In summary, after separating the spin part of the ef- 
fective interaction induced by the transverse gluons, one 
finds that scalar and pseudo-scalar mesons are described 
by the 't Hooft - Gross - Neveu interaction [Eq. ( |2.8D ], 
while vector and axial vect or mesons are described by 't 
Hooft's equation [Eq. dJ), but with 7 = 0]. f\ This is 
one of the principal results of this paper. 

Quantum numbers are associated with these states by 
considering vacuum to meson matrix elements of oper- 
ators with definite quantum numbers (such as iptp and 
tpij^ip). The resulting quantum numbers are consistent 
with the quantum numbers that one can intuitively guess 
by merely considering the degeneracy of states and their 
naive quark model parity properties. Using these rules, 
one obtains the assignments of quantum numbers shown 
in Table I. Of course, while the above assignment of quan- 
tum numbers may appear reasonable, it is still to some 
extent arbitrary, since the model has no rotational sym- 
metry. One should always keep this fact in mind when 



comparing the model results to the experimentally mea- 
sured meson spectra. 

Note that, in the limit of infinite mass for the Aj_ field 
considered here, the model supports only states with he- 
licity wp to h — 1. Tensor and higher spin mesons are 
absent in this limit of the model. 



helicity 


hel. symmetry 


orbital symmetry 







odd 


odd 


0+ 





odd 


even 


0" 





even 


odd 


1+ 


1 


even 


odd 


1+ 





even 


even 


1" 


1 


even 


even 


1" 



TABLE I. Assignment of quantum numbers to the 
states of the model. One of the remarkable results of this 
model is that, despite the rotationally non- invariant formu- 
lation, states with h — and those with h — ±1 are ex- 
actly degenerate (this is true for both vector and axial vector 
mesons) . 



^However, the effective masses of the quarks in the vector 
and axial vector channels are the same as the ones generated 
by the Gross-Neveu interactions in the scalar and pseudo- 
scalar channels. 
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III. NUMERICAL RESULTS 
A. Meson Spectrum 

In the case of vector and axial vector mesons, the nu- 
merical calculations amount to just solving 't Hooft's 
equation, which has been well studied and nothing fur- 
ther needs to be said about this part. For the scalar and 
pseudo-scalar mesons, which are described by a combi- 
nation of the 't Hooft and Gross-Neveu models, the cal- 
culational methods to obtain spectra and wave functions 
are described in the Appendix ^. In both cases (i.e. with 
and without Gross-Neveu interaction) the numerical cal- 
culations were based on variational calculations with a 
basis of Jacobi polynomials, such that the correct end- 
point behavior of the wave functions was manifestly built 
in. typically, a basis size of 5-10 states was already suf- 
ficient to give the lowest eigenvalues with at least four 
significant digits. 

The model contains three free parameters: the 
QCDi+i coupling, the Gross-Neveu coupling and the 
effective (constituent) quark mass. These parameters 
were fixed by using as input the masses of the tt, p 
and p(1450) mesons. The masses of both the p and the 
p(1450) do not depend on the Gross-Neveu coupling, and 
they thus determine the values of the bare parameters in 
the QCDi+i interactions, yielding mq = 288Mey and 
y/G^/iT = 21SMeV respectively. [] The value of the 
renormalized Gross-Neveu coupling is then obtained by 
performinga subtraction at the tt pole as explained in 
Appendix y|. 

The physical string tension can be read off directly 
from the static QQ potential in QCDi+i 



V{x) 



G2 



(3.1) 



For the above model parameters, one thus obtains a phys- 
ical string tension of a = (349Mey)^, which is slightly 
lower than the preferred value used by the lattice com- 
munity [aiatt « [mMeVf]. 

Having determined all free parameters, we can now 
turn to calculate masses of other hadrons. When com- 
paring with experimental data, the comparison is always 
done with isospin one states. The reason for this choice 
is that the model makes a large Nc approximation, and 
it is generally expected that such an approximation is 
better for isospin one states than for isospin zero states, 
since isospin zero states in QCD are affected by mixing 
with gluonic components and such mixing is turned off 
in the large Nc limit. Nevertheless, identifying the model 



spectrum with isospin 1 states is still to some extent an 
arbitrary choice. 

Numerical results for the lightest mesons and their first 
radial excitations are shown in Table II. Given the sim- 
plicity of the model, the agreement with the experimental 
data 1 14 1 is spectacular! 

One may argue which of the mesons in Table II should 
be used as input parameters, particularly since some of 
these states are not believed to be simple quark model 
states but rather complicated "molecule" type states and 
one would not expect that such mesons are adequately 
described by a large Nc model. However, since all 
mesons in Table II are fitted quite well, picking a dif- 
ferent set of input parameters has only a minor impact 
on the quality of the fit. 

However, and this is more important, as we discussed 
in the previous session, there is the issue of assignment 
of quantum numbers, which reflects itself here in pick- 
ing states from the particle data table |0 and com- 
paring the model results to this data. Because of the 
truncation of transverse momenta, the model has less 
degrees of freedom than QCD (or a 3-1-1 dimensional 
quark model). Therefore, not all mesons with spin 
or 1 have a model equivalent. By making (somewhat ar- 
bitrary) assignments of quantum numbers, one does in 
essence pick a subset of mesons from the particle data 
table and fits these mesons. As an example, the way we 
assigned quantum numbers, there are oq and bi mesons 
in the model but no ai mesons. In principle, one might 
thus be tempted to generate a fit to a different sub-set 
of mesons. However, since the above choice of quantum 
numbers seemed to be by far the most reasonable choice, 
no such results will be presented here. Nevertheless, one 
could always keep this in mind as an option. 



meson 




model (MeV) 


exp.(MeV) 


TT 


0" 


139 (input) 


139 


P 


1" 


769 (input) 


769 ± 1 


ao 


0+ 


978 


984 ± 1 


bi 


1+ 


1249 


1231 ± 10 


7r(1300) 


0" 


1339 


1300 ± 100 


p(1450) 


1~ 


1465 (input) 


1465 ± 25 


a*, 


0+ 


1606 ± 2 


7 


bl 


1+ 


1707 


? 



''Note that this quark mass includes already the mass renor- 
malization due to the Gross-Neveu interactions and should 
thus be regarded as a constituent quark mass. 



TABLE II. Calculated masses of the lightest mesons and 
their first radial excitations for each quantum number. Un- 
certainties are indicated only when they exceeded 1 MeV. 

Motivated by the unexpected success of the model for 
calculations of the meson spectrum, we will in the fol- 
lowing consider other observables as well. Of course, be- 
cause of the severe approximations used in constructing 
the model, we will limit ourselves to observables where 
one would expect that the truncation of transverse mo- 
menta is not so critical. Furthermore, we will focus on 
those observables where the LF formalism is helpful in 
simplifying the calculations. 
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B. Meson Wave Functions 

One quantity of great interest, which is easily acces- 
sible in our model, is the so called (twist-2) light-cone 
wave function of mesons, defined by a correlator along a 
light-like direction, i.e. for example 



(Pxdx 



(3.2) 



in the case of the tt, and correspondingly for the p. In 
gauges other than the LF-gauge = 0, a gauge string 
in Eq. (|3.2| ) is understood. LF time x+ is the same for 
the two field operators. Since the definition of these LF 
wave functions involves only longitudinal correlations, it 
seems sensible to consider them in our coUinear model. 



TV i.z 


^\ P 


—LH 

0.6 
0.4 
0^ 


— ,^1 O.Z 3 0.4 5 O.f) 7 OS Q'J 

LO 

as 

\ 0.0 
-0.5 

ao 


0.1 0.2 0.3 0.4 0.5 O.G 0.7 0.8 0.9 

bi ^"^'''^^^-.^ 



0.3 0.4 05 0.6 



X X 

FIG. 2. Light-cone wave functions for the 4 lightest mesons. 

Numerical results for the LF wave functions of the four 
lightest mesons are shown in Fig. ^. The wave func- 
tions for the TT and p are normalized to f^^ and fp re- 
spectively, since this is common in discussions of these 
mesons. No such physically motivated normalization ex- 
ists for the ao and 61 mesons and they are thus normalized 
to dxip'^lx) = 1. The p and bi are just solutions to 't 
Hooft's equation. The boundary behavior (x — > 0, 1) is 
^ ^ x^ [(1 — x)^], where 7r/3cot7r/3 = 1 — m'^n/G^ and 
the number of nodes in the wave function increases with 
the excitation energy of the corresponding meson. 

Obviously, this a crude model and one should be care- 
ful when comparing these results with other calculations 
or with experimental data. In particular, since modes 
with large transverse momenta have been omitted, one 
should only compare to calculations or data at a low mo- 
mentum scale. With these caveats in mind. Fig. ^ shows 
several interesting features, which might be relevant for 
the real world. First of all, the pion wave function is 
nearly flat and it does not vanish at the end-points. This 



is not just a numerical coincidence, as one can see by 
studying vacuum to meson matrix elements of the color 
singlet axial current 



fnP' - m\n). 



(3.3) 



In the chiral limit, d^^jl^ — and thus /„ must vanish 
— unless = 0. Only massless mesons are allowed to 
have /„ 7^ in the chiral limit. In terms of the LF wave 
functions. 



fn 



dx(t>n{x) 



and thus (regardless of the quark mass) 



(3.4) 



(3.5) 



If the TT is the only massless meson in the chiral limit (it 
is in this model), then Eq. ( |3.5D implies 



dx(f>j^{x) — 1, 



(3.6) 



which, together with the normalization condition 
dxtf)"^ {x) = 1 and the Cauchy-Schwarz inequality, im- 
plies 4>-n{x) = 1, i.e. a constant tt wave function. It 
should be emphasized that this argument does not work 
in more than 1+1 dimensions, where one also has to inte- 
grate over fcx- However, in any l-fl dimensional model 
with chiral symmetry and with only one massless meson 
state the wave function of this state must be flat. This 
statement holds even in the presence of higher Fock com- 
ponents, in which case one can also prove that (again if 
there is only one massless meson state) its wave function 
must be purely valence and flat. 

The p meson wave function is peaked in the middle 
and it does vanish at the end points. The main reasons 
for the difference between the tt and p meson wave func- 
tions are the stronger binding and the attractive short 
range interaction, which acts only in the pion channel. 
In a non-relativistic Schrodinger picture, the main differ- 
ence between the tt and the p is that the tt wave function 
gets completely "sucked" into the region of the attrac- 
tive short range interaction, thus leading to very high 
momentum components its wave function. End-points 
of the wave function on the LF, correspond to high mo- 
menta in an equal time picture, which explains why the 
TT wave function is much larger than the p wave function 
near the kinematical end-points (x = and x — 1). 

Scalar and axial vector mesons show similar features 
as the TT and p respectively, except that there is a node 
in the wave function, which arises since these mesons 
are simply the second lowest states in infinite ladders of 
states with an increasing number of nodes — very much 
similar to non-relativistic states in one space dimension. 
The reason for the difference between the end-point be- 
havior of the ao and 61 wave functions is the same as for 
the TT and p. 
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Note that parton distribution functions in this model 
are simply the squares of the distribution amplitudes in 
this model, both because higher Fock components are 
negligible and because there are no integrals over trans- 
verse momenta. However, before comparing the resulting 
parton distributions with experimentally measured struc- 
ture functions, one should bear in mind that the coUinear 
approximation, i.e. neglecting quanta with k± ^ 0, 
makes sense only at low scales. Since it was not clear 
how to determine the starting scale for QCD evolution 
within the model, no attempt was made to compare the 
resulting structure functions to the available experimen- 
tal data. 



C. Meson Form Factors 

In the case of the p meson, the 3 diagrams depicted in 
Fig. H contribute to the elastic form factor. As a side re- 
mark, it should be emphasized that in the LF framework, 
one can usually calculate form factors directly by taking 
matrix elements between the states and one does not need 
to worry about various diagrams (at least for the "good" 
component of the current). However, in the large Nc 
limit, the bound state equations are greatly simplified 
by leaving out all components of the wave function that 
are 1/Nc suppressed. On the LF this not only eliminates 
non-planar diagrams, but also diagrams which are planar 
but which contains fermion lines that go forth and back 
in time. When one calculates form factors, there are dia- 
grams where the photon couples to such backward going 
lines, which become suddenly allowed, even on the LF, 
as long as the -I- component of the momentum transfer 
is nonzero. Examples of such diagrams are shown in Fig. 
|b) and c). 



The expression for the elastic form factor of the p 
meson is identical to the corresponding expression for 
mesons in QCD2 jl^. The reason is that the Gross- 
Neveu interaction neither contributes to the equation for 
the p nor does it contribute to the interaction in the 
photon (vector current!) channel. The form factor for 
mesons in QCD2 has been derived in Rcf. 

F+ = {2p+ - q+) F{q^) = 2p+(i - x) [/""p + r + r] 

(3.7) 



(3.8) 



where x = /p'^ , 

I dz(t>{x + {\-x)z)(l){z) 
Jo 

is the "impulse term" (Fig. ^), 



[x{l -u) + (l- x)zY 

(3.9) 

corresponds to a "vertex-correction" term (Fig. ^d) , and 

r^x^G^[\u[\z t>i^-)'f>i^)Giu;e) ^ (3,^0) 

[x(l — u) -f (1 — x)z] 







is the "exchange current" (Fig. He). G{u;q^) is the inte- 
grated Green's function in QCD2 



G{u-q^ 



G{u,v;q')^J2 



dvG{u, v; (f') 

(t)n{u)(j)n{v) 



(3.11) 



a) 



b) 



(l-x)(l-z) 



/(l-x)z 




c) 



xu 




The fraction x — q^ jp'^ of momentum transfer is deter- 
mined by energy conservation 



2 _ 9 , M 



X \ — X 



(3.12) 



In the case of the tt form factor, the same diagrams con- 
tribute, but there are two additional diagrams which con- 
tain backward moving lines that connect to the meson via 
a Gross-Neveu interaction (Fig. ^ 



FIG. 3. Time ordered diagrams contributing to the vector 
form factor of the p meson: a) impulse term, b) vertex cor- 
rection, c) exchange current. The grey blobs represents the 
meson wave function and the dressed quark photon vertex. 
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a) 



b) 





FIG. 4. Additional time ordered diagrams contributing to 
the vector form factor of the tt. The grey blobs represents 
the meson wave function as the dressed quark photon vertex, 
and the black blob the coupling of a quark to the meson via 
the Gross-Neveu interaction. The slashed fermion line corre- 
sponds to an instantaneous fermion exchange interaction. 



The contribution of these diagrams to [Eq. (3.7)] 
the matrix element of is given by 



FigAa = 2p+x'^Z J du 



M„ 



xu 



(f){xu)G{u] 
(3.13) 



and 



r 1 

FigAb = -2p+x^Z'^ / du 

Jo l~ xu 



Giu;q'), (3.14) 



where Z is a (dimensionful) constant related to the nor- 
malization of the vertex. Its value can be determined 
from 



{x) = -Z I dyGix,y,ti'){^ 



(3.15) 



where is normalized according to jj^ dx(f)'^{x) = 1. Note 
that we distinguished here between the masses of the 
quark and anti-quark, which is not necessary for the tt, 
but which will be useful when we investigate A'-mesons 
later. 

In general, the wave functions and the in tegrals that 
determine the form factor [Eqs. ( |3.8| )-( 3T^ )] have to be 
obtained numerically. However, in the chiral limit the 
pion wave function is constant <j)T,{x) = 1. This yields 
pmp _ ^ Furthermore, f'"'^ and f^"^ exactly cancel each 
other. Finally, with Mq = Mq = M, we find Z = M 
and there is also a partial cancelation between the two 
diagrams in Fig. ^ yielding 

/•I ][,p 

FigA = -2p+x du G{u]q^) 

Jo 1 - M 



■ p'^x 



duG{u; q^) 



(3.16) 



where we used that the solutions to 't Hooft's equa- 
tion satisfy (equal mass case) du(f>n{u) / u{l — u) — 
jl^ Jp du(j)n{u). Adding up the various pieces and using 



furthermore [from Eq. ( 3.12| )] that (7+ = p+ (i.e. a; = 1) 



for vanishing tt mass, we thus find the remarkably simple 
result for the elastic vector form factor of the pion 

F^{q^) = l-q^ f duG{u-q^)^Y.4^^ (3-17) 

where gvi"^) = /q^ du(l>n{u). In particular, we find for the 
radius 



rms-ra 



(3.18) 



i.e. for the parameters used above one obtains a numer- 
ical value of 



0.605/m, 



(3.19) 



which is only slightly smaller than the rms-radius for 
the p meson in the same model \/{r^)^ ~ 0.625. The 
explanation for this near equality of rms radii is that 
both form factors are dominated by the rho meson pole 
(note that sf^jnip = 0.627/m). [| The experimental 
value for the rms-radius of the pion is only slightly larger 
\J {r'^)-K,exp = 0.663 ± 0.006/m. However, one should not 
be too impressed by this excellent agreement, since the 
p mass is an input parameter and since the rms radius 
of the pion is dominated by the p meson pole. The in- 
teresting aspect of the model in this respect is that it 
reproduces vector meson dominance. 

While the form factors of the tt and p mesons are rather 
similar for low momentum transfers, there is a qualita- 
tive difference between them at large momentum trans- 
fers. The elastic form factor of the p meson follows the 
asymptotic behavior of form factors in QCD2, which is 
given by 



1 



(Q^) 



1+/3 



(Q2 ^ 00), 



(3.20) 



where (3 is the same exponent that also appears in the end 
point behavior of the wave function. For the parameters 
used in the present fit {ml k, G^/tt) one finds (3 « 0.5 . 



*Note that we use here the 3+1 dimensional relation between 
the slope of the form factor and the rms radius since this is 
supposed to be a model for 3-f 1 dimensional QCD. 

''Of course, for the p meson it would be the lj meson pole, but 
in the large Nc limit the p and omega mesons are degenerate. 
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Naively, from Eq. ( 3.171 ) one would expect F^^^Q'^) ~ 
as predicted by naive power counting. However, 
^„ /i^gp(n)(7\/(n) diverges logarithmically and a more 
careful analysis yields 



,lnQ2 



(Q' -> oo). 



(3.21) 



This logarithmic growth of Q^FTr{Q^) arises because we 
considered the limit of an infinitely heavy gluon in the 
collin ear QCD model. One should thus not take Eq. 
(3.21) literally at very large Q^, but only in an inter- 
mediate momentum range, where it makes sense to in- 
troduce an infinitely heavy gluon as an effective degree 
of freedom. 




0.05 0.1 0.15 0.2 0.25 
^2 




Q Q 

FIG. 5. Numerical results for the pion form factor [Eq. 
(3.17)]. The experimental data is from Ref. [15]. 

While the details of the form factor analysis at large 
are model dependent, there are again some qualitative 
features which might be important also in full QCD^+i. 
Since there is an attractive short range interaction acting 
only in the tt channel, the wave function of the tt (in an 
equal time picture) contains much more high momentum 
components than the p. As a result, the elastic form 
factor of the tt should fall off much slower than the form 
factor of the p. 



a flavor independent coupling of the quarks to the trans- 
verse component of the gluons and if one strictly takes 
the limit of infinite mass for the transverse gluons. The 
results of a fit, where TOs as the only new parameter is 
shown in Table III. Using the K mass as an input, one 
finds nis ~ 502Mey, which is quite heavy compared to 
a light quark mass of ruq = 278MeF (the fitted value 
from the non-strange sector). As a result, the masses of 
the K* and the (p turn out to be much too high. Alter- 
natively, if one uses the K* mass as as input parameter, 
1^ the K turns out much too light. 

However, one should keep in mind that one should con- 
sider the Gross-Neveu interaction of the model an effec- 
tive interaction, which is obtained by integrating out high 
energy degrees of freedom. In real QCD, these degrees of 
freedom do not have infinite mass and one would there- 
fore expect some scale (and thus also some fiavor-) depen- 
dence of the effective interaction. Since we are not able 
to perform this elimination procedure explicitly, we must 
also regard the Gross-Neveu coupling for strange quarks 
as an independent parameter in this phenomenological 
model. Corresponding results are shown in the forth col- 
umn of Table HI (fit 2). The strange quark mass in this 
case (from fitting the mass of the K^) turns out to be 
TOs = 408Mey, which is much more reasonable com- 
pared to niq than the above fit. The K mass is fitted by 
construction and yields the subtraction constant for the 
Gross-Neveu interaction, ^ which we used to calculate 
the mass of the lightest strange meson with 0'*' quan- 
tum numbers. Compared with the lightest meson with 
these quantum numbers that can be found in the particle 
data table, the model calculation gives a much too small 
mass. This is to some extent surprising, since the model 
did so well for the non-strange 0^ meson. However, we 
will not attempt to explain this discrepancy here, since 
the quark model of these meson is to some extent still 
very speculative |ll6|]. 



meson 




fit 1 


fit 2 


exp.(MeV) 




0^ 


494 ±5 


494 


494 


K* 


1" 


986 


892 


892 


<i> 


1" 


1198 


1012 


1020 


if* (1430) 


0+ 


1180 ± 10 


1129 


1429 ± 6 



D. Strange Quarks 

The extension of the model to include strange quarks 
is straightforward and any details will be omitted here. 
The only issue are the numerical values of the model 
parameters, i.e. which of the parameters one should take 
from the non-strange fit. 

One might argue that the only new parameter is the 
strange quark mass, since the gauge coupling is obviously 
flavor independent and since the Gross-Neveu coupling 
also turns out to be flavor independent if one starts from 



^°Both the K* and the (\> do not depend on the Gross- 
Neveu coupling and are thus better suited for determining 
the strange quark mass independently from the Gross-Neveu 
coupling. 

^^In fit 1 in Table III, the numerical results for the K and the 
/s'*(1430) have large errors associated with them, since those 
calculations do not subtract at the i('-pole and one therefore 
has to deal numerically with divergent quantities. 
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TABLE III. Comparison of numerically calculated meson 
masses to experimental data. In fit 1, all parameters other 
than the strange quark mass are taken from the non-strange 
fit. The Kaon mass is taken as an input parameter. In fit 
2, the Gross-Neveu coupling is also allowed to differ from the 
coupling in the non-strange sector. 



The light-cone wave functions for mesons containing 
strange quarks are quite similar to wave functions of non- 
strange mesons and there is thus no reason to display 
them here. The only difference is that the wave-functions 
are no longer symmetric under x —^ 1 — x, but they are 
slightly shifted such that the s quark tends to carry more 
momentum than the light quark, as one would naively 
expect. 

With the high quality electron beam now available at 
Jefferson Laboratory, there will be attempts to measure 
the electro-magnetic form factor of the if-meson over a 
wide momentum range [ p^ . 

In order to calculate this observable from our coUinear 
QCD model we concentrate on fit 2 (Gross-Neveu cou- 
pling adjusted), since it is crucial to have both the correct 
ii'-meson mass as well as the (/)-meson mass in order to 
obtain a reasonable model for the Kaon form factor. Ap- 
plication of above formulas for the form factor is tedious 
but straightforward and the result of such a calculation 
is shown in Fig. M. 




FIG. 6. Numerical results for the kaon form factor. Dashed 
line: strange form factor, dotted line: up form factor, full line: 
electro-magnetic form factor. 



A comparison with Fig. || shows that the form factors 
of IT and K mesons are rather similar. At first, this seems 
surprising, but the physics of this result becomes clear 
when one decomposes the form factor of the K'^ into the 
contributions from the u and the s quark. The various 
form factors now find a simple interpretations in a quark 
model: Clearly, the u form factor falls off faster than the 
s form factor, since the heavier s quark is more localized 
than the u quark. Furthermore, when comparing the u 
form factor in a if with the u form factor in a tt we find 
that the u quark is more localized in the tt than in the 
K. There are two reasons for this trend. First of all, 
the Gross-Neveu coupling is weaker in the K and thus 
there is a little less attraction in the K compared to the 
TT (this compensates perhaps for the increased reduced 
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mass). Furthermore, while the mean separation between 
quark and anti-quark changes perhaps only little when 
going from the tt to the K, the center of mass gets shifted 
towards the s quark, which means that the u quark gets 
"pushed out" a little. The combination of a slightly faster 
(compared to the tt) falloff of the u form factor and the 
slower falloff of the strange form factor (compared to the 
d in the tt) then results in a net electro-magnetic form 
factor for the K meson which is very similar to the tt 

form factor. 

For the rms radii one finds Vvl)i^ ~ 0.490 ± .005 fm 
and \/Jr^}K ~ 0.665±.005/m resulting in a net (electro- 
magnetic) rms radius of -y/ {r^)K ~ 0.612 ± .005/m (ex- 
periment: 

x/{r^)K+ ~ 0.58±.04/m [|l8|), which is almost 
identical to the rms radius of the tt •^/(r^)^ w 0.605/to 
in this model. For the neutral Kq, the model predicts 
{f'^)K° ~= 0.067 ± .005 fm^, which is consistent with the 
experimental result (r2)^o « -0.054 ± .026 fm^ @. 

While numerical model results are consistent with 
available experimental data ||l^, the error bars of the 
latter are large and more precise data, over a larger mo- 
mentum range, would be very useful |jl7| . 

IV. SUMMARY 

We have considered the collinear model for QCD where 
all modes with nonzero transverse momentum (relative to 
an arbitrary, but fixed, direction) are neglected. Such a 
model has not been derived from QCD, and should thus 
be regarded as a purely phenomenological model. Within 
this model, we focussed on the limit where the effective 
mass of the transverse gluon becomes very large, but with 
an appropriately rescaled coupling to the quarks. As a 
result, one obtains a model that resembles the 't Hooft 
model but with helicity degrees of freedom for the quarks 
and point-like helicity-dependent quark-antiquark inter- 
actions that resemble the chiral Gross-Neveu model. This 
limiting case of collinear QCD is solvable for large Nc. 

Given the ad-hoc truncation of the degrees of freedom, 
one would be tempted to dismiss such a model as not 
being a very useful model for QCD^+i. However, as the 
physical observables presented in this paper indicate, the 
model yields a remarkably good description for many ex- 
perimental data. 

First of all, the spectrum of vector and axial-vector 
mesons turns out to be rotational invariant (i.e. helicity 
independent !). 

Secondly, with only three free parameters available 
[which are fixed to reproduce the tt, p and yo(1450) 
masses] one obtains an extraordinarily good result for 
the masses of other non-strange mesons. One of the most 
surprising result in the calculations of the spectrum was 
the excellent agreement with the mass of the ao(980), 
which is often difficult in quark models and which is be- 
lieved to possess a large KK component (which is not 



there at large Nc). We have no explanation for this sur- 
prisingly good agreement of the model results with the 
experimental data. This agreement might be accidental 
and it does also depend on which physical mesons are 
identified with meson states within the model. Since the 
model lacks transverse momenta, the model contains of 
course much less states than for example the quark model 
and therefore one must pick a subset of states with which 
the model states are identified. While the choice made 
in the paper seemed by far the most reasonable, there is 
some uncertainty in the assignment of quantum numbers 
since the model lacks full rotational invariance. 

The spectra of mesons containing strange quarks, were 
less well reproduced. If one uses the same coupling in the 
strange and non-strange sectors, then it turns out to be 
impossible to fit the K and K* masses at the same time. 
While the fit improves considerably, when one allows the 
Gross-Nevcu coupling for strange quarks to differ from 
the coupling for light quarks, the overall fit for mesons 
with strange quarks never reaches the quality of the fit 
for non-strange mesons. Specifically, the = 0"*" meson 
is not very well reproduced, which might be a hint that 
the agreement in the non-strange sector is to some extent 
accidental. 

Motivated by the success of the meson spectrum calcu- 
lations, we considered other observables as well. An ob- 
servable which is particularly easily accessible in the LF 
framework are the LF wave functions of hadrons. The 
wave function for the tt turned out to be almost com- 
pletely flat [(t)-^{x) = 1 in the chiral limit]. In contrast, the 
wave function of the p turned out to be strongly peaked 
in the middle, which is more reminiscent of a weakly 
bound state. Physically, this difference did arise because 
the model has an attractive zero-range interaction which 
acts only in the tt channel. A similar picture arises when 
one considers the LF wave functions for the ao and bi 
mesons. Both wave functions have nodes in the middle, 
which reflects the parity of these mesons, but only the 
wave function for the ao is non-vanishing at the bound- 
ary. 

One feature which makes this model particularly use- 
ful for form factor calculations is that it is fully covariant 
under boosts in the longitudinal direction, which makes 
the extraction of form- factors frame independent. De- 
spite the differences in their LF wave functions, the p 
and TT meson turn out to have almost the same rms ra- 
dius (i.e. slope of the vector form factor. The main 
reason for this result is vector meson dominance for the 
vector form factor at small momentum transfers, which 
is also a feature of collinear QCD. At large momentum 
transfers, the p form factor falls off much faster than the 
TT form factor. Roughly speaking, a flat LF wave function 
corresponds to an equal time wave function which has a 
very large high momentum component (the end points 
of the LF wave function correspond to infinite momenta 
for the constituents in an equal time framework) thus 
resulting in less wave function suppression for the large 
momentum transfer form factor. 
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In order to obtain a reasonable fit for the spectrum of 
mesons containing strange quarks, it turns out to be nec- 
essary to re-fit the Gross-Neveu coupling. After having 
done this one obtains an electro-magnetic form factor for 
the K meson, which is very similar to the form factor of 
the TT. This result is "explained" by s form factor which 
falls of slower and a u form factor that falls of faster than 
d and u form factors in the tt so that the net difference 
between K and tt nearly vanishes. The predicted charge 
radii are consistent with experimental data for K'^ and 
mesons. 



(a) 



(b) 
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APPENDIX A: ON THE RENORMALIZATION 
OF THE VERTEX MASS 

As long as the same cutoffs are used in a covariant cal- 
culation and in a LF calculation |^ the current quark 
mass of the covariant calculation and the vertex mass of 
the LF calculation usually are the same. Chiral symme- 
try breaking manifests itself through the renormalization 
of the kinetic mass term in the LF calculation p9|-|2^ . 

The situation changes when the LF calculation is done 
using a Tamm-Dancoff approximation. This fact is best 
illustrated by considering a concrete example, such as 
a fermion interacting with some boson field. For sim- 
plicity, let us assume that the fermion self-energy has 
only a scalar piece and let us furthermore assume that 
this scalar piece is momentum independent ^ — in other 
words, let us assume that 



= (5m, 



(Al) 



where 5m is a constant. 

Without truncations of the Fock space, the three point 
interaction for the physical fermion (=solution to the 
Hamiltonian) gets modified through Feynman diagrams 
that include self-energies connected to the vertex by an 
instantaneous fermion line (Fig. |^). 



^■^Under these circumstances, the only difference between a 
covariant calculation and a LF calculation is that zero-modes 
are omitted in the latter. 

^^An explicit model with such features is provided by a 
model where the boson mass goes to infinity and where the 
coupling is rescaled such that a nontrivial mass renormaliza- 
tion remains. 



(C) 



FIG. 7. Diagrams contributing to the three-point vertex for 
a physical fermion. (a) bare vertex, (b) self-energy insertion 
on the outgoing fermion line, (c) self energy insertion on the 
incoming fermion line. 

Note that under the above assumptions about the self- 
energy S, there are no further diagrams contributing to 
the physical vertex. 

For the case of a (transverse) vector vertex the bare 
vertex yields a matrix element for helicity fiip transitions 



r 



flip _ 
hare 



1 



P 



(A2) 



where = p+ — p'^ is the momentum transfer at the 
vertex. The two diagrams with a self energy insertion in 
the incoming and outgoing external lines yield, respec- 
tively 



"J- out 



1 



q+ 

dm 1 



(A3) 



so that 



flip 
bare 

V ^ 



' in 

dm ( 1 



flip 



(A4) 



This example shows that, if one imposes a Tamm-Dancoff 
(TD) truncation such that the 5m correction to the ver- 
tex can no longer be generated through higher order Fock 



Other cases are analogous. 
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components, then one needs to renormalize the vertex 
mass accordingly in order to obtain the same physical 
results as without the TD truncation. Hence, with a TD 
approximation in place, the numerical value of the ver- 
tex mass need no longer agree with the numerical value 
of the current mass. In particular with a TD approxi- 
mation, there is no reason that the vertex mass should 
vanish in the chiral limit. This is in contradistinction to 
LF calculations without TD approximations [pO[. 



APPENDIX B: THE 'T HOOFT 
NEVEU MODEL 



GROSS 



Because the Gross-Neveu interaction is quadratic in 
the non-dynamical component of the fermion field V'-j 
the constraint equation becomes non-linear and a canon- 
ical LF formulation of models with such an interaction 
is somewhere between difficult and impossible. Usually, 
one has to use some indirect methods for solution. In this 
Appendix, we will use covariant methods, supplemented 
by the LF solution to QCZ3i+i, to construct a solution 
to the Gross-Neveu model with additional QCDi+i in- 
teractions. 

Let us consider the 1+1 dimensional model 



£ = V - ^tr [F^^F^ 

2_ 

27V 



(Bl) 



where ip has color components (fundamental represen- 



tation), Ffj^i, = d^A^ - dtyAf^ + i-^ [A^, A^] is the gauge 
field and the limit — s- cx) is implied. Obviously, for 
7 = one obtains the 't Hooft model (QCDi+i in the 



limit A^ ^ oo) |10|, while for G = the model reduces to 
the chiral Gross- Ncvcu model |l^ . Because of the four 
fermion interaction, £ has to be renormalized. ^ 

Before we der ive the bound state equation correspond- 
ing to Eq. (Bl), let us illustrate the techniques in the 
simpler case G — (the chiral Gross-Neveu model): First 
of all, the scalar tadpole gives rise to a "constituent mass" 
M for the fermion. Secondly, since we are considering the 
N oo limit, the meson spectrum can be determined by 
summing up a (geometric) series of bubbles, yielding for 
the pseudo-scalar two-point function 



1 



tr 



1 



(2^2 - (g^_M) (^-y-Af ) 



(B2) 



^^We closely follow Ref. |12|, where more details can be 
found. 

^''Up to an overall normalization factor, which is not relevant 
for determining bound state masses. 



The integral in Eq. (B2) diverges logarithmically. We 
renormalize by subtracting at = (strictly speaking 
we first introduce a regulator, then subtract and then 
send the regulator to infinity), yielding 



dx 



Jren ^TT Jq IVP - p'^x{l - x) 



(B3) 



where the renormalized (finite!) coupling is related to the 
bare coupling (zero!) via (cut-offs not explicitly shown) 



1 

7 
1 



(2^2 " [4-M) 



tr- 



1 



dx 



7 2ti Jq x{1 — x) 



(B4) 



The pseudo-scalar bound state (in the following referred 
to as the pion) mass is obtained in terms of the renor- 
malized coupling be searching the pole of D^ip^). In 
practice, one often turns the argument around and uses 
the pion mass as an input parameter, which yields the 
renormalized coupling as a function of the pion mass 



= 



Ir 



AP - fj,lx{l - x) ' 



(B5) 



Since the 't Hooft model is most conveniently solved in 
the LF framework, we will also formulate the 't Hooft 
- Gross - Neveu model using LF quantization. As a 
warmup exercise, we thus reconsider the Gross-Neveu 
model in this framework as well. Using standard canoni- 
cal quantization plus renormalization of the kinetic mass 
plus (as explained in Appendix ^) renormalization of the 
vertex mass one finds the (still ill defined!) bound state 
equation for pseudo-scalar mesons 



APij{x) jAP f dy ip{y) 
x{l-x) x{l-x)J 2Try{l-y) 



(B6) 



(for scalar mesons, where the wave-function is odd under 
X 1 — X, another term contributes which is omitted 
here for simplicity). Note that we set the vertex mass 
equal to the kinetic mass in Eq. (B6). As is explained 
in Appendix ^ the vertex mass can acquire nontrivial 
renormalization, if certain corrections are suppressed by 
the approximations used (here: fermion tadpoles con- 
nected to the bare vertex by an instantaneous fermion 
line). The fact that we set it here equal to the kinetic 
mass term is arbitrary and irrelevant, since the interac- 
tion term is also multiplied by the (still undetermined) 
bare coupling constant. But, in order to achieve maxi- 
mum similarity with the covariant approach, we set the 
vertex mass equal to the kinetic mass at this point. 

Eq. ( p6| ) yields a divergent solution, since il^ix) ^^^^ 
const, and thus J dy'ip{y)/y{l — y) diverges. This diver- 
gence is exactly the same divergence that one obtains in 
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the covariant approach as well (before coupling c onst ant 
renormalizatfon) as one can see by reducing Eq. (B2) to 
a parameter integral, yielding 



dx 



7 2tt Jo x{l-x)M'^-p'^x{l~x) 



(B7) 



and the whole divergence can be avoided provided one 
renormalizes properly. Here we will use a trick |^ to arrive 
at a renormalizcd bound state equation: the operator 
identity 



(B8) 



where m is the current quark mass, together with Lorentz 
covariance, implies for vacuum to meson (on shell) matrix 
elements 



p~(0|?/'7+75V'|n,p) ^ ■m{0\-ipj5il;\n,p), 
i.e. in terms of the wave-function 

M 



(B9) 



^J■ 



m I dx — - — 



x) 



M^) (BIO) 



The reason M and not m appears in the LF expression 
for (0|'!/'i75^/'|7r(j))) is again the Tamm-Dancoff expansion 
used. We can use this identity to express the divergent 
integral in Eq. ( |B6|) by a convergent integral, yielding the 
renormalizcd (but non-hermitian, since the Hamiltonian 
depends on the eigenvalue) bound state equation for the 
pion channel 



fi^ipix) 



x{l — x) mx{l — x) 



j'i*ivl (Bl 



1) 



This "miracle" is possible since the current quark mass 
itself is cutoff dependent and goes to zero (for fixed phys- 
ical masses) as the cutoff is sent to infinity. In fact, it 
goes to zero in such a way that the ratio 7/771 stays finite, 
which motivates the definition 



7^^ 

/ ren 



(B12) 



A consistency check of this result is obtained by explicitly 
solving the renormalizcd integral equation. From Eq. 
(Bll) one finds 



ip{x) 



trenr" 



— x{l — x)fl'^ 

Consistency requires that 



dy 
2tt 



(B13) 



'^''A similar trick has been used by T. Heinzl (talk given 
at the "1997 workshop on LF quantization" in Les Houches, 
France). 



LF 2 



dx 



27r - x{l- x)^l'^' 



(B14) 



which we recognize as the pole condition for Dt^{p^) [Eq. 
(p3|)] with 7ren replaced by 7^^. 

Let us now proceed to the full 't Hooft ~ Gross -Neveu 
model (Bl). A calculation analogous to the covariant 
procedure discussed above can again be obtained by sum- 
ming up a chain of pseudo-scalar bubbles, but now each 
bubbles is dressed by the QCDi+i interactions as well as 
by the Gross - Neveu (tadpole) self-energy. |^ 




FIG. 8. In the pseudo-scalar two-point function, the qq 
propagation within each bubble is governed by the QCDi+i 
Kernel in planar approximation (-^ 't Hooft model). 

In practice, these bubbles can be evaluated by inserting 
a complete set of ('t Hooft-) meson states 



7 27r ^ ^ 



(B15) 



where the /i„ and gp{n) are the masses and pseudo-scalar 
couplings of the n-th meson solution for the 't Hooft equa- 
tion. The masses are most conveniently calculated in 
terms of the LF wave-functions 



gp{n) — M / dx 



1 



1 - X 



tpnix). 



(B16) 



The sum in Eq. ( B15D diverges, but since QCDi^i is 
superrenormalizable, this is only the free field divergence. 
We renormalize again by performing a zero momentum 
subtraction, yielding 



1 



1 P^gpjnf 



(B17) 



where 



^^The fact that the fermion lines within the Gross - Neveu 
tadpoles are dressed by QCDi+i interactions only changes 
the self-mass by a finite amount. Since the induced mass is 
only a bare parameter anyway, this additional renormalization 
of a parameter that is already renormalized has no further 
consequences. 
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1 ^ I gpjnf 

Iren 7 „ 



(B18) 



is finite. For the scalar channel we are not free to renor- 
malize independently, since chiral symmetry demands 
that the scalar and pseudo-scalar couplings are the same. 
This yields 

n-i/ 2n 1 1 sr 9s{nf_ 



Iren 27r ^ fll [jll - p^Y 



(B19) 



where 



gs{n)^Mj^ dx (^^ - Mx) (B20) 
are the scalar couplings and 

1-K ^ hi 



1 1 



accounts for the (finite) differences in the zero momentum 
subtraction in the scalar and pseudo-scalar channels. For 
practical calculations, the following representation is also 
useful 



Eqs. (B17) and (B19) or (B22) are perfectly suitable 



for a numerical determination of the spectrum of the 't 
Hooft - Gross - Neveu model. However, since we want 
to determine LF wave-functions and parton distribution 
functions, we also need to understand how to renormal- 
ize the LF wave equation for this model. The canonical 
procedure (modulo renormalization of both kinetic and 
vertex masses as explained above) yields the (still ill de- 
fined) LF bound state equation for this model 



dy- 



1 dy 



■>P{y) 



x{\ — a;) 

a;(l-a;)yo 27r y(l - y) 
7M2(1-2x) dy{l^2y)ij(y) 



{x- 



x{l — x) 



2tt y{l-y) 



(B23) 



In the pseudo-scalar channel, we can use the same pro- 
cedure that we introduced in connection with the Gross 
- Neveu model, namely replacing the pseudo-scalar cou- 
pling by the vector coupling, which yields the renormal- 
ized LF bound state equation for the 'Hooft - Gross - 
Neveu model in the pseudo-scalar channel (i.e. only even 
wave functions) 



lYipix) 



x{l - x) Jo 
x{l - x) Ja 2tt 



V'(y) 



(x - y)2 
ijiy). (B24) 



It is quite ea sy to verify consistency of the renormalized 
LF equation ( B24 ) with the more conventionally obtained 
result (B17). Using the Green's function for the 't Hooft 
equation. 



Gix,y,p')^J2 



1pn{x)'lpn{y) 



(B25) 



one can rewrite Eq. (B24) in the form 

^(x)= / dyG(x,y,/.2)^ff^ / -^(z). (B26) 
Jg y(l - V) Jo 27r 

Upon integrating over x and using the fact that the so- 
lutions to 't Hooft 's equation satisfy 

f dxMx) = f dx4r^ (B27) 
Jo Jo ^ X) 

provided we identify with 



we recover Eq. (B17) 



^ren • 

For scalar mesons, we have not been able to derive a 
renormalized LF bound state equation. However, this 
does not prevent us from calculating the LF wave func- 
tions of sca lar r esonances, using the following trick: ^ 
From Eq. (B23) one notes that the wave functions for 
scalar states satisfy 



tpnix) cx / dyG{x,y,fil) 



I'M 



M 



1 



(B28) 



where is the invariant mass of the bound state. Thus 
one can first obtain the bound state mass from the poles 
of D„{p^) [i.e. from Eq. (B22)] and then one determines 
Tpnix) (up to a normalization fac tor) by solving the linear 
equation corresponding to Eq. (B28), i.e. one solves 



(Mo 



M ( - 

X 



1 



(B29) 



for ip{x), where Hth is the 't Hooft Hamiltonian. 
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